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FOREWORD 


This  report  was  prepared  by  Technology  Incorporated,  Dayton,  Ohio, 
in  compliance  with  United  States  Air  Force  Contract  Number  AF  33(657)- 
8073,  initiated  by  the  Retardation  and  Recovery  Branch,  Flight  Accessories 
Laboratory,  Directorate  of  Aeromechanics,  Aeronautical  Systems  Division. 
The  original  Air  Force  Project  Engineer  was  Mr.  Clint  Eckstrom;  in  the 
later  phases  of  the  program,  Mr.  Lawrence  L.  Watson  assumed  this  position. 

Technology  Incorporated  inaugurated  the  program  on  15  February  1962 
and  completed  its  investigation  on  15  November  1962  \uider  the  supervision 
of  Mr.  William  B.  Walcott,  Project  Engineer.  Other  staff  members  con¬ 
tributing  significantly  to  the  project  were  Mr.  J.  Patrick  Ray,  who  performed 
the  statistical  analysis,  and  Mr.  Krishan  Joshi,  who  conducted  the  literature 
survey. 


ABSTRACT 


A  study  was  conducted  to  determine  the  effects  of  changing  scale  upon 
drag  coefficient,  filling  time,  peak  opening  force,  and  stability  for  single, 
unreefed  textile  parachute  canopies.  The  investigation  was  confined  to  Flat 
Circular,  Extended  Skirt,  Ringslot,  Ribless  Guide  Surface,  Circular  Flat 
Ribbon,  and  Conical  Ribbon  parachutes  operating  in  the  subsonic  flow  regime 
at  altitudes  below  20,  000  feet.  As  an  integral  part  of  the  study,  a  survey  of 
the  existing  literature  and  test  data,  a  compilation  of  all  pertinent  data,  and 
recommendations  for  future  experimental  investigations  were  made.  Scaling 
equations  with  associated  95-percent  confidence  intervals  were  developed  for 
the  drag  coefficient  of  the  Flat  Circular  and  Extended  Skirt  parachutes  through 
the  use  of  a  multiple  regression  analysis  which  indicates  the  significant  varia¬ 
bles  and  their  functional  forms.  Because  the  available  data  was  poorly  dis¬ 
tributed,  the  equations  will  have  to  be  used  circumspectly  to  avoid  misleading 
conclusions . 


Publication  of  this  technical  documentary  report  does  not  constitute  Air  Force 
approval  of  the  report's  findings  or  conclusions.  It  is  published  only  for  the 
exchange  and  stimulation  of  ideas. 
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SECTION  I 
INTRODUCTION 


Suitable  scaling  equations  for  textile  parachutes  would  substantially 
reduce  the  testing  necessary  to  define  their  operating  characteristics.  In¬ 
stead  of  using  current  methods  of  testing  and  modifying  full-scale  parachutes 
until  the  desired  results  are  obtained,  scale  models  in  conjunction  with  the 
scaling  equations  would  suffice;  then,  if  desired,  a  full-scale  model  could  be 
employed  to  verify  the  results. 

This  study  was  intended  primarily  to  determine  whether  sufficient  data 
exists  which  could  define  appropriate  scaling  equations  for  the  selected  para¬ 
chute  types.  Data  voids  and  future  test  needs  were  to  be  noted.  If  sufficient 
information  were  available,  the  scaling  equations  were  to  be  derived. 

To  determine  the  scaling  equations,  it  was  decided  at  the  outset  to  use 
a  multiple  regression  analysis,  a  statistical  technique  believed  to  be  novel 
in  its  application  to  the  parachute  field.  Several  reasons  dictated  the  use  of 
a  regression  analysis  rather  than  the  more  conventional  engineering  methods 
commonly  applied  to  a  study  of  this  nature.  Any  other  type  of  analysis  could 
not  handle  the  many  variables  involved.  These  variables  were  numerous 
because  the  data  for  determining  the  scaling  equations  were  collected  from 
many  different  tests;  for  example,  if  sufficient  data  were  available  over  the 
desired  canopy  area  range,  other  variables,  such  as  different  canopy  loadings, 
and  different  chute  porosities,  would  also  be  present.  Furthermore,  use  of 
a  regression  analysis  would  enable  determination  of  the  significance  of  a  given 
variable,  the  functional  form  of  the  variable,  the  goodness  of  fit  of  the  equa¬ 
tion,  and  the  confidence  limits  for  the  final  results. 


SECTION  II 

PARACHUTE  SCALE  EFFECTS 


A  .  Basic  Concepts  of  Scale  Effects 

A  parachute  scale  effect  can  be  defined  as  any  change  in  a  parachute's 
operating  or  performance  characteristics  caused  solely  by  a  change  in  canopy 
size.  Such  a  change  is  brought  about  because  change  in  canopy  size  can  alter 
certain  geometric  and  aerodynamic  properties  of  the  parachute.  These  effects 
can  be  significant  in  some  situations  or  virtually  negligible  in  others,  but  are 
always  present  to'  some  degree  in  any  scale  model  testing.  As  in  any  scale 

Manuscript  released  by  the  author  September  1962  for  publication  as  an  ASD 
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model  testing,  if  this  effect  cannot  be  assessed,  the  test  results  can  prove  to 
be  useless.  The  main  difficulty  in  parachute  scaling  is  the  impossibility  of 
assessing  this  change  theoretically. 

The  change  in  parachute  performance  characteristics  as  effected  by 
scaling  manifests  itself  as  a  change  in  parachute  drag  coefficient,  stability, 
filling  time,  or  peak  opening  force.  If  these  variables  are  defined  as  the  de¬ 
pendent  variables,  then  the  problem  is  to  define  their  change  in  terms  of 
variations  of  the  independent  variables  caused  by  a  change  in  parachute  size 
(scale).  The  most  obvious  independent  scaling  variable  is  parachute  canopy 
area;  however,  there  are  several  other  independent  variables  which  reflect 
changes  in  size  that  may  be  significant.  Reynolds  number  and  Froude  number 
are  such  independent  variables  since  by  their  definitions  they  vary  with  a 
change  in  size.  Other  independent  variables  which  can  be  effected  by  a  change 
in  size  are  the  number  of  gores  and  the  ratio  of  their  number  to  the  canopy 
diameter.  It  is  obvious  that  one  or  both  of  these  changes  will  occur  with  a 
change  in  parachute  size. 

However,  a  variance  in  gore  number  does  not  necessarily  constitute  a 
scaling  variable.  There  are  three  questions  which  must  be  satisfied  before 
a  variable  can  be  considered  a  scaling  variable.  These  questions  can  be 
enumerated  as  follows: 

1)  Does  the  variable  change  with  a  change  in  scale? 

2)  Does  the  variable  affect  the  performance  of  the  parachute? 

3)  Is  the  only  way  to  assess  the  effect  of  the  variable  to  consider  it 
simultaneously  with  a  change  in  canopy  size  (scale)? 

Due  to  the  manner  in  which  the  number  of  gores  in  a  parachute  vary  with  a 
change  in  canopy  size,  the  answers  to  all  three  questions  must  be  in  the  affirm¬ 
ative,  and  for  question  3  this  is  reasonable  when  it  is  considered  that  a  1-foot 
nominal  diameter  scale  model  of  a  50-foot  nominal  diameter  chute  cannot 
possibly  have  the  same  gore  arrangement  as  the  model.  The  50-foot  nominal 
diameter  chute  will  have  approximately  50  gores  while  for  practical  reasons 
the  1-foot  diameter  model  would  probably  be  constructed  with  10  to  20  gores. 

To  ascertain  the  effect  of  the  change  in  gore  arrangement  it  is  necessary  to 
consider  it  in  conjunction  with  the  change  in  scale  since  the  scale  model  will 
not  be  constructed  with  either  the  same  number  of  gores  or  the  same  ratio 
of  number  of  gores  to  nominal  diameter  as  the  full  scale  chute.  At  the  same 
time  the  effect  of  this  change  in  gore  arrangement  cannot  be  assessed  inde¬ 
pendently  since  scale  model  chutes  and  large  diameter  chutes  will  probably 
never  be  constructed  with  the  same  gore  arrangement  due  to  the  practicalities 
of  construction. 

Other  variables  which  affect  parachute  performance  but  are  not  scaling 
variables  since  they  do  not  satisfy  question  3 — since  their  effects  can  be 
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assessed  independently — are; 

a)  porosity 

b)  suspension  line  length/ nominal  diameter 

c)  number  of  vertical  ribbons 

d)  number  of  horizontal  ribbons 

e)  space  between  vertical  tapes/gore  width  at  base 

f)  open  spacing /ribbon  width 

g)  canopy  loading 

h)  riser  length/nominal  diameter 

i)  vent  area/nominal  surface  area 

j)  velocity 

k)  dynamic  pressure 

Even  though  these  above  listed  variables  cannot  be  considered  as  scaling 
variables,  they  still  may  be  present  in  the  scaling  equations  due  to  their 
presence  in  an  interaction  term. 


A  special  class  of  Independent  variables,  of  which  Reynolds  number  and 
Froude  number  are  examples,  is  defined  as  an  interaction  term.  If  a  term  is 
differentiated  with  respect  to  the  variable  denoting  canopy  size  and  other  varia¬ 
bles  are  still  present,  the  term  is  defined  as  an  interaction  term.  For  example. 


which  upon  partial  differentiation  with  respect  to 


dRe  _  _£V 
dDo  ~  ^ 


becomes 


From  this  rather  simple  example  it  can  be  seen  that,  if  Reynolds  number  is 
significant  for  scaling,  the  variables  >  ,  ^  ,  and  v  are  also  significant. 

Therefore,  even  though  these  variables  are  held  constant  from  the  full  scale  to 


the  scale  model  tests,  the  value  of  the  variable  will  determine  the  amount  of 
change  in  the  parachute's  operating  characteristics  due  to  scaling.  From  the 
preceding  discussion  it  can  be  concluded  that  almost  any  variable,  for  example, 
velocity,  might  possibly  play  an  important  role  for  scaling  if  it  appears  in  a 
significant  interaction  term.  An  interaction  term  which  was  found  to  be  sig¬ 


nificant  for  the  definition  of  the  scaling  effects  on  filling  time  and  opening 
pDq 

shock  was  -  .  It  is  possible  that  interaction  terms  involving  other  vari¬ 


ables  may  also  be  significant  even  though  none  were  found. 


In  summary  then,  the  aerodynamic  and  geometric  parameters, 

a)  nominal  surface  area, 

b)  number  of  gores, 

c)  number  of  gores /nominal  diameter, 
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d)  Reynulda  number,  and 

e)  Froude  number 

can  be  assumed  to  be  scaling  variables  and  the  other  parameters,  Itemized 
previously,  will  enter  into  the  scaling  equations  only  if  they  are  present  in 
interaction  terms. 

B.  Discussion  of  the  Literature  Surveyed 

Literature  was  surveyed  extensively  to  determine  whetlier  sufficient 
test  data  were  available  to  permit  an  analysis  leading  to  the  definition  of  the 
appropriate  scaling  equations.  The  data  collectetl  during  this  survey  is  sum¬ 
marized  and  presented  in  Appendix  II.  The  survey  was  limited  to  Solid  Flat 
Circular,  Extended  Skirt,  Ringslot,  Ribless  Guide  Surface,  Flat  and  Conical 
Ribbon  parachutes  opi-rating  at  indicated  airspeeds  below  300  knots  and 
altitudes  h?as  than  20,000  feet. 

Since  the  data  were  collected  from  a  great  many  different  reports,  there 
were  many  variables  present  which  could  affect  the  different  performance 
characteristics  of  the  parachutes  and  yet  not  be  true  scaling  variables.  Since 
this  was  the  case,  it  was  necessary  to  consider  several  more  independent 
variables  than  would  normally  be  required  to  arrive  at  the  scaling  equations. 
The  independent  variables  whose  values  were  considered  necessary  were 
S  o  I  1$  /  Do  •^g'^Do,lf/Dg,Sy/So,^,W/So(V,Q,Rd,M,  and  Fr.  For 
this  reason,  a  data  point  would  be  discarded  unless  a  reasonable  and  logical 
assumption  could  be  made  of  the  value  of  the  missing  variable  or  variables. 

A  more  complete  discussion  of  the  literature  survey  is  inchjded  in  Appendix  II 
with  the  data  summary. 

It  was  felt  that  sufficient  data  was  available  to  delennine  the  regression 
equations  lor  drag  coefficient  for  the  Solid  Flat  Circular,  Extended  Skirt, 

Flat  Ribbon,  Ringslot,  and  Ribless  Guide  Surface;  parachutes.  Als.u,  there 
appeared  to  In-  sufficient  data  for  filling  time  anti  peak  opening  force  for  both 
the  Flat  Circidar  and  Extended  .Skirt  parachutes.  However,  for  reasons 
which  will  be  (explained  later,  the  only  attempts  to  derive  scaling  equations 
which  yielded  reasonable  results  w«;re  for  the  Flat  Circular  and  Extended 
Skirt  parachut»;s  for  drag  coefficient. 

OthtM-  than  tin;  data  ftjr  the  above  mentioned  parachutes,  the  available 
data  were  eltln-r  insufficient  <n’  nonexistent.  Data  for  jiny  of  the  performance 
charac  ti;  ristics  for  the  Conic  .il  Ribbon  p.iraidiute  wen-  so  meager  that  any 
attempt  to  summari/.i-  and  present  them  would  serve  no  ustdul  purpose.  Since 
no  common  description  «)r  definition  of  stability  has  been  developed  in  the 
parachute  realm,  tin-  available  data  in  the  liti-rature  could  not  be  utilized. 

For  full -sc;i  1«;  pa  r.ichutes ,  the  period  of  oscillation  or  .'iinount  of  oscillation 
was  acquin-d.  h'oi-  scale;  model  wiinl  tunin-l  t«-sts,  howevt;r,  pitching  moment 
and  the  stabh-  angle  of  attack  were  olit.iiin-d.  Mecause  of  this  .and  the  fact 


-1 


that  much  of  the  stability  data  appeared  to  have  been  observed  rather  than 
measured,  any  attempt  to  define  the  scaling  equations  or  to  present  the  data 
in  a  systematic  form  was  abandoned. 

C.  Utilization  of  Regression  Analysis 

Because  of  the  scatter  of  data  caused  by  the  many  variables  in  the  col¬ 
lected  data,  a  multiple  regression  analysis  was  selected  as  the  technique  most 
capable  of  defining  the  desired  scaling  equations.  Figure  1,  for  example, 
shows  the  wide  scatter  of  drag  coefficient  values  for  various  nominal  surface 
areas  for  the  Circular  Flat  parachute.  This  figure  shows  that  there  are  other 
variables  in  addition  to  canopy  size  which  influence  the  value  of  the  drag  co¬ 
efficient,  for,  without  these  other  variables,  the  scatter  for  any  nominal  sur¬ 
face  area  would  be  appreciably  less.  Therefore,  it  was  necessary  to  try  not 
only  to  extract  the  scaling  effects  from  the  data  but  also  to  eliminate  the 
variations  due  to  nonscaling  variables,  such  as  suspension  line  length  and 
velocity. 


As  an  example,  assume  that  a  regression  equation  has  been  developed 
for  the  prediction  of  drag  coefficient  which  involves  only  two  term3--one  term 
is  a  scaling  variable  (In  Sq)  and  the  other  term  is  not  (In  V)  i.  This  equation 
can  then  be  used  to  extract  from  the  observed  data  points  the  variation  con¬ 
tributed  due  to  velocity  and  thus  allow  the  utilization  of  many  more  data  points 
than  if  all  the  observed  data  points  had  to  be  at  the  same  velocity.  Therefore, 
it  can  be  seen  that  the  regression  equation  defines  the  variation  due  to  scaling 
(In  Sq  )  and  also  makes  possible  the  utilization  of  data  in  which  nonscaling 
variables  affecting  performance  are  also  varied. 


The  multiple  regression  analysis,  the  theory  of  which  is  outlined  in 
Appendix  I,  offers  several  attractive  features  for  an  analysis  of  this  type. 

If  the  data  to  be  analyzed  by  a  regression  technique  are  reasonably  well  dis¬ 
tributed  over  the  ranges  of  interest  for  the  independent  variables,  it  is  pos¬ 
sible  to  determine  the  following: 

a)  how  closely  the  developed  regression  equation  is  describing  the 
variations  for  the  observed  data  points, 

b)  the  significance  of  any  given  independent  variable  by  using  appropriate 
statistical  techniques, 

c)  the  functional  form  of  a  significant  independent  variable,  and 

d)  the  confidence  limits  for  the  predicted  dependent  variable. 

Two  methods  are  generally  employed  in  choosing  the  regression  model 
to  represent  the  functional  relation  between  the  dependent  and  independent 
variables.  One  is  to  plot  the  sample  data  as  2-dimensional  scatter  diagrams 
of  the  dependent  variable  versus  each  of  the  independent  variables  with  the 
other  independent  variables  held  constant  to  give  some  indication  of  what  the 
individual  terms  of  the  model  should  be.  The  second  method  is  to  base  the 
functional  form  of  the  regression  model  on  prior  knowledge  and  analytical 
considerations  of  the  factors  involved. 

The  first  method  was  not  adaptable  for  utilization  in  the  present  study 
since  this  study  is  based  on  data  gathered  from  various  sources  and  collected 
under  radically  different  experimental  conditions.  It  was  impossible  to  find 
enough  data  collected  under  similarly  controlled  experimental  conditions  to 
yield  meaningful  2-dimensional  plots.  The  second  method  was  only  partially 
feasible  since  a  thorough  review  of  the  pertinent  literature  yielded  very  little 
information  applicable  to  the  combined  effects  of  the  independent  variables 
considered  in  this  study. 

Another  source  of  information  indicative  of  the  independent  variables 
and  cross-products  of  independent  variables  (interactions)  which  are  signifi¬ 
cantly  affecting  the  dependent  variable  is  a  statistical  computation  procedure 
known  as  analysis  of  variance,  a  technique  for  analyzing  the  results  of  a 
designed  experiment.  A  designed  experiment  is  one  planned  in  advance  with 
the  levels  and  combinations  of  levels  of  the  independent  variables  chosen  so 
that  the  maximum  desired  information  on  the  factors  may  be  elicited.  The 
analysis  of  variance  method  is  based  on  a  partitioning  of  the  variance  of  the 
observations  on  the  dependent  variable,  each  part  measuring  the  variability 
attributed  to  some  specific  source.  The  method  provides  the  quantities 
necessary  for  testing  the  significance  of  this  variability.  Obviously,  this 
technique  was  not  applicable  in  the  present  study  due  to  the  poorly  distributed 
data  on  the  independent  variables. 
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It  was  decided  to  use  a  multiple  regression  technique  in  building  a  regres¬ 
sion  model  by  a  process  of  trial  and  elimination.  Application  of  this  technique 
to  the  parachute  field  is  believed  to  be  unique.  (Appendix  I  discusses  the 
utilization  of  the  IBM  7090  in  the  development  of  this  model.  )  The  data  for 
this  development  were  gathered  and  processed  in  the  form  of  four  dependent 
variables  and  approximately  12  independent  variables  considered  to  have  the 
greatest  effect  on  the  former.  A  regression  model  for  drag  coefficient  in¬ 
volving  linear  terms  of  the  independent  variables  was  first  fitted  but  rather 
unsuccessfully.  Other  functional  relations  — such  as  models  composed  of  the 
squares,  reciprocals,  or  natural  logarithms  of  the  independent  variables  — 
were  all  fitted  with  varying  degrees  of  success.  One  functional  relation  would 
fit  reasonably  well  for  one  type  of  parachute  but  would  not  fit  at  all  for  some 
other  type  of  parachute. 

With  the  possible  exception  of  differences  in  the  coefficients,  approxi¬ 
mately  the  same  functional  form  for  a  given  dependent  variable  should  fit  all 
different  types  of  parachutes.  For  each  respective  dependent  variable,  the 
approach  employed  to  obtain  similar  functional  relations  for  all  types  of  para¬ 
chutes  was  the  following;  Various  transformations  on  the  independent  varia¬ 
bles —  transformations  such  as  squaring,  computing  the  reciprocal,  squaring 
the  reciprocal,  extracting  the  square  root,  computing  the  natural  logarithm, 
squaring  the  natural  logarithm,  and  computing  the  reciprocal  of  the  natural 
logarithm  of  "1"  plus  the  variable — were  performed;  furthermore,  cross- 
product  terms— such  as  the  number  of  gores  divided  by  the  nominal  diameter 
and  dynamic  pressure  times  nominal  surface  area— were  computed.  Each  of 
these  terms  was  considered  to  be  an  "independent"  variable.  The  correlation 
coefficients  for  all  possible  pairs  of  these  "independent"  variables  were  cal¬ 
culated,  and  the  highly  correlated  terms  were  removed.  Iterative  multiple 
regression  analyses  were  run  on  the  remaining  terms  until  all  those  not  sta¬ 
tistically  significant  had  been  deleted.  The  resulting  functional  relations 
are  presented  in  Tables  1  and  2. 

Table  1  presents  the  drag  coefficient  prediction  equations  for  the  five 
types  of  parachutes  which  had  sufficient  data  to  attempt  a  regression  analysis. 
Performing  the  transformations  on  the  independent  variables  presented  in  the 
table  and  then  combining  all  terms  results  in  a  predicted  drag  coefficient. 

The  degree  of  fit  (R^)  of  the  tabulated  (observed)  values  in  Appendix  II  for 
the  prediction  equations  is  indicated  in  the  parentheses.  The  same  procedure 
follows  for  Table  2  which  presents  peak  opening  shock  and  filling  time  pre¬ 
diction  equations  for  the  Circular  Flat  and  Extended  Skirt  parachutes.  For 
a  given  dependent  variable,  dissimilarities  in  these  functional  relations  for 
the  various  types  of  parachutes  must  be  attributed  to  differences  in  the  avail¬ 
ability  of  data  for  the  various  parachutes. 

In  an  effort  to  explain  as  much  of  the  data  variation  as  possible,  several 
interaction  terms  which  appear  in  various  theoretical  parachute  analyses  were 
tried  in  the  regression  models.  Except  for  Reynolds  number  and  Froude 
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T.ible  1 


Prediction  Equations  for  Drag  Coefficient 


Ribless  Guide  Surface 
(r2  =  .  9253) 

Elxtended  Skirt 
(R2  =  .8833) 

Circular  Flat 
(R2  =  .8635) 

R  in  slot 
(R2  =  .839  ,) 

Ribbon 
(R-  =  .6351) 

2. 7367 

8. 5736 

14.  80  *0 

-I6l.0:«b 

17.  1465 

7895  In  Sp 

f.  2676  IftSg 

13l6in  So 

-.0760  Sq 

t  3.  7779  In  Sg 

t.0271Sp 

r.  7067  In 

-H.  2727y?7“ 

-.  1024  W/S, 

-  0553  Sp 

-.0251 

f.  3985  F^ 

♦  5 .  boss  tn  F^ 

♦  .  379oin  w/Sg 

-  1.  8017^57 

-.5919  l/lW/Sp) 

-.1312  1/F, 

♦  .  7340  F^ 

•■.0007i/(w/So)* 

t9  8516  1  /Sp 

f.O'in  i/(w/Spi* 

-.  2195Mp/D„ 

*1.  1263  (InF,)* 

7492  in  V 

0305  l/tW/Sp) 

-.0008  W/Sp 

r.  8292  In  W/Sp 

-.6929  l/(W/S^ 

r. 0015  V 

4858  In  V 

ctl.  6853  t/V 

5ti00  InV 

..  1852  (In  W/S,)* 

rb. 1*86  Ny 

-  4265  In  Rp 

f.0022  V 

-2.  6to8  InV 

t.  0831  l/(W/Sj,)* 

tU83.7375  l/N, 

1601  k 

-.0285  1/9* 

-1.  1010  In  V 

52.5409  l/k 

-25.  5995  (10^)  1/flp 

-.0556  9 

12.2819  In  X 

t82.6927  (10*0)  l/H* 

t. 2306  1/9 

+5.6517  yr 

..3235l/(l«/Op)* 

-.0245  1/9* 

+.0089  I/M 

-.8069  (10-5)  I/M* 

Prediction  Equations  for  Filling  Time  and  Opening  Shock 


1  Filling  Time 

Opening  Shock  | 

Circular  Flat 

Extended  Skirt 

Circular  Flat 

Elxtended  Skirt 

(r2  =  .5121) 

(R^  =  .4137) 

(R'i  ••  .  5564) 

(r2  .  .8695) 

-  .  8279 

-5oo4. 0513 

210.  7387  (I02) 

-1213.  9166  (10^) 

+  3142.  5799  ^Oq/v’ 

3M0  1  /0 

-4757.9144  (10*)  , 

+9771.  1968  (10)  InN, 

-  1881.  1376  (10^)  ,* 

+  1573.  1123  (10^)  l/H# 

+  9780.  2371  (10*)  /Dp/V  * 

+  2861.  5800  I/M 

* . 0029  V 

>474  2279  In  R« 

-366.7599  (10^)  i/lnda-q) 

-2520.  9976  (10)  InS, 

- . 0010  \/0 

t467. 0649  In^ 

-7.  3409  Sp 

+  1081.  3l'35  (105)I/n, 

t . 004  1  k 

-.  3352  (10-*)  l/n* 

+  3.9308  (10-5)  R, 

.51.8082  I/M* 

♦  303.  8672  Sy/Sp 

-5664. 4910  (10*5) |/R,» 

-a. 1269  k 

+  1238.  1923  (10)  InM 

+  .  3037  (lO-O)  1/.* 

-1694.  6732  (10^)  n 

-947.  1344  Nj 

+  35.  3660  1  /  X 

-  3889  (10-5) 

-.  7334  9 

.1 1  36.  0969  (109)  ,* 

+  17.4472  Inq 

+  .0015  q* 

■3783  4468  (I0)n0,/V  * 

+3  8789  Ff 

+  lb  5306  In  Sp 

•  114.  3176  If  Ft 

+4.  7390  Sg 

-4028.  1048  (10)  N,/0a 
-4051.  9106  (10*)  l/Sp* 
-138.  1699  I/IW/S,)* 
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number,  however,  the  only  interaction  term  found  to  be  significant  was  _ 

V* 

which  was  for  filling  time.  This  does  not  necessarily  indicate  that  these  other 
variables  may  not  be  significant;  for  this  could  have  resulted  from  the  ob¬ 
served  data  used  in  the  regression  analyses.  It  is  of  interest  to  note  the  b 

A  Df\ 

value  (coefficient)  for  - for  the  Circular  Flat  parachute  listed  in 

V 


Table  2.  This  value  with  its  associated  95-percent  confidence  limits  will 

g 

encompass  the  constant, — (/•q=.002378)  .  which  is  given  in  the  "USAF 
Parachute  Handbook"  (reference  30). 


It  was  initially  anticipated  that  the  same  variables  with  similar  functional 
forms  and  possibly  just  slightly  different  coefficients  would  be  significant  for 
all  the  parachute  types.  That  the  regression  models  did  not  substantiate  this 
anticipation  was  apparently  due  to  the  poor  distribution  of  the  available  data. 
Nevertheless,  it  is  believed  that  adequate  data  will  prove  the  hypothesis.  At^ 

least  it  should  be  noted  that  the  scaling  variables  (  In  S^,  Fr  ,  In  Fr ,  Fr ,  In  Frj  ) 
for  the  Circular  Flat  parachute  drag  coefficient  are  similar  to  the  scaling 
variables  (In  S<, ,  In  Fr,  Fr,  l/Fr)  for  the  Extended  Skirt  parachute  drag  coefficient. 
Furthermore,  these  two  regression  models  were  obtained  from  larger  amounts 
of  data  than  were  available  for  the  other  types  of  parachutes. 


Scaling  Equation  Theory 


If  the  prediction  equations  as  presented  in  Tables  1  and  2  are  used, 
the  scaling  equations  can  be  derived  as  follows;  The  only  variables  of  in¬ 
terest  in  these  prediction  equations  are  the  scaling  variables;  therefore,  for 
practical  purposes  the  other  variables  are  disregarded.  The  reasoning  that  ■ 
substantiates  this  statement  is  that  a  scale  effect  is  interpreted  as  a  change 
in  the  dependent  variable  (drag  coefficient,  etc.  )  due  to  a  change  in  only  the 
independent  scaling  variables.  As  a  change  denotes  an  increment,  the  non¬ 
scaling  variables  can  be  eliminated  since  they  don't  vary  through  this  interval. 
Assuming  the  prediction  equation  to  have  five  terms  of  which  the  first  two 
and  the  last  are  scaling  variables  (b,  and  bg  X^j  )  , 

=«0  *  +‘■4*4, + 

then,  the  incremental  change  in  the  dependent  variable,  as  predicted  by  the 
regression  equation,  is 

Yg-Y,  =  b,  (  X,2  -  X„)  +  bg  (  Xgg-Xg,)  +  bgCXgg-Xg,)  =  b,  K,  +  bg  Kg  +  bg  Kg 

Thus,  (  Yg-Yj )  is  the  incremental  change  in  drag  coefficient,  filling  time,  or 
opening  shock  caused  by  a  change  in  scale  and  this  derived  equation  is  then 
the  scaling  equation. 


With  the  change  in  the  dependent  variable  known,  it  is  also  desirable  to 
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determine  the  confidence  limits  for  this  change.  To  predict  confidence  limits, 
the  variance  of  the  predicted  change  must  be  estimated. 


To  estimate  the  variance  of  the  difference  between  two  predicted  mean 
values  of  the  dependent  variable  when  one  or  more  of  the  independent  varia¬ 
bles  remains  constant,  let  it  again  be  assumed  for  purposes  of  illustration 
that  the  fitted  regression  equation  has  five  terms: 


9  = 


/V  A 

The  estimate  of  the  variance  of  (  —  Y^ )  is  desired  when,  for  example, 

and  remain  constant,  i.e.,  X31  ~  ^32  ^4/ “^43* 


Now, 

/\  /V  r  /\  /\  2 1  r  /V  A  T 2 

Var  (Y2-Y,  )=e[(Y2-Y,)  J-[e(Y2-Y,)J  , 

^2  ~  ~  ^  ^I2~  ^  ^  ^2^  ^21^  ^5  ^  ^52  ~  ^  ^2*^2^  ^5 


Thus , 


Var  (Y2-Y,)  =  E[(b,K,  +  b2K2  +  b5K5)^]-[E  (  b,  K,  +  b^  +  bg  K^)  ] 

Considering  these  two  terms  separately, 

E  [{ b,  K,  +  bg  Kg  +  bg  K5)^]=E  [  bfKf+  bg^Kg^^  b/Kg^  2  b,  bg  K,  Kg  +  2  b,  b^  K,  K5+  2 bg  b^  Kg  ] 

+  2K|Kg[oj2+^l^+ 2  K|K5j^0j5+/3|^5j+  2  K^K^  [^25  +  ^g^g] 


and 


[E(b,K,+b2Kg  +  b5K5)]  =[k,^,+K2^2+K5^5]  =  K,  +  Kg  ^g +K5 

+  2  K,  K2;3,^2+  2  K,K5  /35+  2  KgK^/Sg  ^5 

Taking  the  difference  of  these  two  expressions  yields 

Var  ( Yg  -  Yj )  =  +  Kg^or^^ +  Kg  2KjKgCr|g+ 2K|Kg  o-^g  +  2 K^Kg  o-g^  , 

where  o-j^=Var(b.)  and  j ,  =  Cov  { bj  ) 
would  be  cstiiTiated  by 


s"9_9  =  s^"[k 

2  I 


, " C„  +  Kg"  Cgg  +  K/C55  +  2  K,  Kg  C,g  +  2  K,  Kg  C,5  +  2  KgKg  Cgg  ] 


wliere  c  .  and  C.  are  defined,  as  before,  to  be  elements  of  the  inverse 
1 1  I) 

matrix  C  . 
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The  predicted  change  due  to  scaling  with  associated  95-percent  confidence 
limits  is  then, 

1^2 -^,'±196  9^2-?, 

E,  Limitations  of  Regression  Models  for  Development  of  Scaling  Equations 

There  are  several  reasons  why  the  regression  equations,  other  than  the 
equations  for  the  Circular  Flat  and  Extended  Skirt  parachutes,  were  not  used 
to  develop  scaling  equations.  Presented  in  Table  3  are  the  ranges  of  the  in¬ 
dependent  variables  of  the  observed  data.  If  the  regression  equations  are 
used  for  prediction  purposes  for  values  of  the  independent  variables  outside 
these  ranges,  a  correct  result  would  be  purely  accidental.  However,  not  only 
is  it  necessary  to  consider  the  range,  but  also  the  distribution  of  the  observed 
data  with  respect  to  the  independent  variables  in  this  range.  Of  primary 
importance  for  this  investigation  is  the  distribution  with  respect  to  canopy 
nominal  surface  area.  These  distributions  for  surface  area  are  presented 
in  Figures  2  through  8  and  are  typical  of  the  distributions  for  the  other  in¬ 
dependent  variables.  As  can  be  seen,  the  ranges  for  the  independent  variable 
presented  in  Table  3  can  be  quite  misleading  since  only  extreme  values  are 
indicated.  It  is,  therefore,  necessary  to  consider  the  distribution  of  ob¬ 
served  data  points  for  every  scaling  variable  before  the  scaling  eqxiations 
can  be  developed  and  used. 


Table  3 

Ranges  of  Independent  Variables 


F«rftebut« 

Dapendent 

So 

1.^0 

iSj/D. 

s  /s 

>1 

W/So 

V 

q 

B.  I  10*® 

D 

Pr 

Typ* 

Varlablea 

(rt.^) 

•t.F'l/ZlIjO 

(rtVrt’-.in) 

(P.f) 

(ft/aec 

li»r) 

B 

0.1<« 

0.616 

1.00 

0.0 

0.0 

0 

o.^ee 

10. >. 

0  .1383 

.0094 

.229 

Flet 

to 

17,671 .6 

to 

1.810 

to 

20.70 

ID 

1.3 

to 

0.25 

to 

426 

to 

51.60 

to 

2jj.6 

m 

to 

28.650 

to 

.2445 

►-■c* 

•  o 

Circ  ular 

^b2.0 

0.700 

B 

Bl 

m 

60 

0  .256 

168.'=’ 

M 

5.45 

to 

to 

IH 

to 

to 

to 

to 

to 

to 

iiiim 

e!>«.o 

0.850 

B 

B 

B 

260 

0.499 

467.6 

204.4 

.4236 

15.58 

1.91C 

0.600 

0.5007 

c.o 

0.0 

10 

0.119 

9.8 

0.108 

0.129 

0.0088 

,374 

Cd 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

Extended 

5319.73 

1.400 

17.96 

0.1 

0.1 

275 

10.33 

111.7 

12.55 

10.053 

0.100 

15.75 

SKlrt 

■•f 

715.6 

0.699 

C.500 

0.0 

ra 

30 

mi 

150.0 

24.9 

30.56 

0.1340 

4,274 

and 

to 

to 

to 

to 

to 

to 

Bll 

to 

to 

to 

to 

to 

5230.0 

1.000 

0.9376 

0.09 

o.uoe 

135 

0.736 

504.1 

194.8 

280. 0 

0.4681 

15.05 

3.66 

1 .000 

1.345 

0 

c 

(■.078 

9.9 

O.lOB 

0.747 

0.0087 

0.511 

RlbtKjn 

Cl 

to 

to 

to 

to 

tTj 

to 

to 

Lt 

to 

to 

to 

to 

129.2 

1  .400 

7.412 

1 .430 

0.010 

JO.O 

70.53 

514.0 

295  .C 

35.40 

0.4627 

30.93 

46  .0 

1.0 

.3346 

s.o 

O.IO 

9.8 

0.111 

C.727 

.0088 

0.503 

Ring  Slot 

No 

to 

to 

to 

Value 

to 

to 

to 

to 

to 

to 

to 

920.6 

1.535 

1.964 

0.0(W 

B8B 

46  .20 

584.0 

364.0 

25  .55 

0.510 

36.63 

0.61*. 

B 

i  .2  TO 

B 

B 

0.161 

0.469 

>”^.009 

0.629 

Gulae  Sorfase 

Cl, 

to 

tc 

to 

It 

Bl 

to 

to 

to 

CRltiess ) 

48. ]» 

13.5t 

2  75 

39.90 

0.471 

45.0 

•Note.  .‘L  rsTtge  •*  ■jeua'etrlc  ^.roi-lty  ^'Iven  In  per  cetiL. 
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Extended  Skirt  Filling  Time  and 
Opening  Shock  Data  Distribution 


Figure  3 

Circular  Flat  Filling  Time  and  Opening 
Shock  Data  Distribution 


Figure  6 

Ringslot  Drag  Coefficient  Data 
Distribution 


toot* 


NOMINAL  SURFACE  AREA  (1**1 

Figure  4 

Extended  Skirt  Drag  Coefficient  Data 
Distribution 


Figure  7 

Ribbon  Drag  Coefficient  Data 
Distribution 


Figure  8 

Ribless  Guide  Surface  Drag 
Coefficient  Data  Distribution 
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As  an  obvious  example  of  why  a  scaling  equation  was  not  developed, 
virtually  all  of  the  collected  data  for  filling  time  and  opening  shock  tor  the 
Circular  Flat  parachute  fell  in  one  area  range.  Therefore,  any  predictions 
involving  canopy  areas  outside  this  range  would  be  unwarranted  extrapola¬ 
tions.  It  should  also  be  noted  with  reference  to  Figure  4  that  data  for  the 
Extended  Skirt  parachute  drag  coefficient  was  nil  for  the  smaller  nominal 
surface  areas.  For  this  reason,  the  restriction  on  the  use  of  the  developed 
scaling  equation  in  this  canopy  area  range  should  be  carefully  observed. 
Also,  reference  to  Appendix  II  will  show  that  the  distribution  of  the  observed 
data  for  Ng  /  Dg  for  the  Extended  Skirt  is  for  practical  purposes  much 
narrower  than  indicated  in  Table  3. 


The  moat  subtle  and,  consequently,  most  difficult  reason  to  recognize, 
which  prevents  the  development  of  scaling  equations  from  most  of  the  regres¬ 
sion  models,  is  the  warped  or  unreasonable  surface  shape  predicted  by  the 
regression  models.  In  other  words,  it  is  very  possible  for  the  regression 
model  to  have  a  very  high  multiple  correlation  ( R  )  ,  but  due  to  the  poor  dis¬ 
tribution  of  the  observed  data,  be  capable  of  predicting  accurately  only  for 
the  observed  data  points. 

9000 

1000 


traa  »  130  POINTS 
■—  •  55  POINTS 
■  >19  POINTS 

'■)  «  I  pf)i  N  r 
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20 
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RFYNOLOS  NUMBER 

Figure  9.  Ribbon  Parachute  Observed  Data  Distribution  for  Scaling  Variables 


To  clarify  this  point,  the  regression  model  for  the  Circular  Flat  Ribbon 
parachute  will  be  used  as  an  <!xample.  Figure  9  presents  the  observed  data 
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distribution  for  the  Ribbon  parachute  as  a  function  of  the  scaling  variables, 
canopy  surface  area,  and  Reynolds  number.  These  variables  were  chosen 
to  demonstrate  the  data  distribution  since  both  were  significant  scaling  varia¬ 
bles  as  predicted  by  the  Ribbon  parachute  regression  model.  Figure  9  shows 
that  the  observed  data  is  very  highly  concentrated  in  three  or  four  regions. 
For  comparison  purposes.  Figure  10  presents  the  observed  data  distribution 
for  the  Circular  Flat  parachute  for  its  significant  scaling  variables,  canopy 
surface  area,  and  Froude  number.  The  Circular  Flat  parachute  data  can  be 
seen  to  be  much  more  evenly  distributed  as  a  function  of  its  scaling  variables 
than  the  Ribbon  parachute  data  is. 


- - -  •  (20  POINTS 
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■  •  10  POINTS 
A  •  S  POINTS 
o  ■  I  P0l|(T 
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Figure  10.  Circular  Flat  Observed  Data  Distribution  for  Scaling  Variables 


To  demonstrate  the  effect  that  the  difference  in  data  distribution  can 
have,  assume  that  the  drag  coefficient  is  predicted  for  both  types  of  para¬ 
chutes  for  a  75-square  foot  surface  area,  10^  Reynolds  number,  and  1. 35 
Froude  number  (Reynolds  number  and  Froude  number  are  compatible  so  that 
the  same  condition  is  being  investigated  for  both  parachute  types).  It  will  be 
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found  that  the  regression  model  for  the  Circular  Flat  parachute  will  predict 
very  well  and  the  regression  model  for  the  Ribbon  parachute  will  predict  very 
poorly.  It  can  be  seen  by  reference  to  Figures  9  and  10,  however,  that  the 
point  at  which  the  drag  coefficient  prediction  is  being  attempted  falls  among 
clusters  of  observed  data  for  both  parachute  types.  To  demonstrate  the 
cause  of  this  apparent  paradox,  the  Ribbon  parachute  regression  model  can 
be  differentiated  with  respect  to  canopy  surface  area,  which  yields, 


3S^ 


.0553  + 


3.7799 


1.8017  _  9  6516 


Evaluation  of  this  differential  through  a  canopy  surface  area  range  of  40  feet 
to  120  feet^  yields  the  slopes  for  the  predicted  surface  as  presented  in  Figure 
11.  Since  the  observed  data  for  the  Ribbon  parachute  is  very  heavily  con¬ 
centrated  around  the  canopy  surface  areas  of  55  feet^  and  110  feet^,  it  can 
be  surmised,  with  reference  to  Figure  11,  that  the  prediction  surface  is 
passing  through  the  observed  data  at  these  points  and  dropping  down  to  an 
inflection  point  at  a  canopy  surface  of  75  feet^.  This  in  actuality  is  what 
is  occurring  and  this  is  demonstrated  in  the  3-dimensional  sketch  of  the  pre¬ 
dicted  surface  presented  in  Figure  12.  The  observed  data  points  are  rep¬ 
resented  by  the  dots  and  the  prediction  by  the  cross  in  this  figure.  It  can  be 
seen  from  the  sketch  of  this  surface,  how  a  regression  equation  can  fit  the 
observed  'data  very  well  (high  )and  still  not  be  appropriate  for  prediction 
purposes  at  any  other  points.  This  is  probably  the  most  forceful  argument 
that  can  be  made  for  the  necessity  of  a  well  designed  experiment  which  results 
in  a  favorable  distribution  of  observed  data  points. 


Figure  11.  Slope  of  Drag  Coefficient  Surface  for  Ribbon  Parachute 
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Figure  12.  Drag  Coefficient  Surface  Predicted  by  Regression  Equation 


F.  Scaling  Eqviations  for  Circular  Flat  and  Extended  Skirt  Parachutes  Drag 
Coefficient 


This  subsection  presents  the  information  necessary  to  derive  the  scaling 
effects  and  associated  confidence  intervals  on  drag  coefficient  for  the  Circular 
Flat  and  EIxtended  Skirt  parachutes  and  gives  an  example  using  the  Circvilar 
Flat  parachute.  The  reasons  for  selecting  only  these  two  regression  models 
to  develop  scaling  equations  are  as  explained  previously. 


Presented  in  Table  4  are  the  scaling  variables  and  coefficients  for  the 

Circular  Flat  and  Extended  Skirt  parachutes  which  may  vary  for  a  change  in 

scale  if  the  other  independent  variables  are  held  constant.  Table  5  presents 

the  associated  c..  and  c..  values. 

II  ij 


Table  4 


Table  5 


Scaling  Variables  and  Coefficients 


Flat  Circular 

Extended  Skirf**  | 

Variable 

bj 

Variable 

In  So 

0.  1316 

mfm 

0.2676 

yPT 

-11.2727 

0.  7067 

InF. 

5.6058 

■GIH 

0. 3985 

F,  - 

0.7340 

KM 

-0. 1312 

(lnF,r 

1.  1263 

-0.  2195 

^Information  to  be  uaed  only  in  the 
nominal  aurface  canopy  area  range 
of  100  feet^  to  5300  feet^. 


Inverse  Matrix,  C,  Values 


Flat  Circular 


Extended  Skirt 


C„  = 

Ci2  = 
ic 

iCj3  - 

'»4  = 
1C, 5  = 
22  = 
C23  = 
IC24  = 


0.E257 
-2.8230 
2. 1354 
0. 1666 
0.  3342 
148.9704 
-76.4048 
-10.0956 


C25  =  -14.9632 


Cj3  - 

C34  = 
C35  = 
C44  = 
C45  = 
C55  = 


1.  9722 
5.0599 
7.9017 
0. 7060 
0.9812 
1.5620 


'"11 

C12 

Cl  3 
Cl4 
Cl5 
C22 
C23 
C24 


1.6231  C 


25 


1.6305 


6.9017  C33  =  1.0708 

-0.2483  C34  =  -0.9737 
0.1911  €33  e. 0.6395 
0.1528  €44=  1.1988 

32.6942  €45=  0.5544 

-2.  7609  C55  a  0.  3892 
2.7680 


S^p  =  0.  006508 


S^E  =  0-  003025 
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To  predict  the  change  in  drag  coefficient  for  the  Circular  Flat  parachute 
for  a  change  in  size  from  a  nominal  diameter  of  two  feet  to  a  nominal  diameter 
of  fifty  feet,  assuming  velocity  eqvial  to  32.  2  feet/second,  the  procedure  is  as 
follows: 

First,  with  Dq  =  50  and  Dq=  2  ,  the  variables  listed  in  Table  4  would  yield 
the  following  values: 


0 

50 

2 

In  So 

7.  58248 

1.14473 

A" 

.89586 

2.00312 

In  Ff. 

-.21995 

1.38941 

.80256 

4.01249 

(In 

.04838 

1.93046 

where 


Now, 


K|-  I  X,2“  X,|) 

where 

and 

Xgg"  ^21^ 

where 

and 

*^3~*  ^32”  X31) 
where 
and 


X, 2=7. 58248 
X„  =  1.14473 


X22=.  89  5  86 
X2,=  2.00312 

X32  =-21995 
X3,  =  1.38941 


K4  -  (  X^2  ^41  * 
where 
and 

*^5~*  ^52“  ^51  ^ 
where 

and 


X^2=  80256 
X^,  =  4.01249 

Xg2=  04838 
Xg,  =  1.93046 


ACp  =  0I3I6(7,58248-II4473)-II.2727(089586-2.00  3I2) 

+  5.6058  (-  0.21995  -  1.38941  )  +  0.7340  (0.80256  -  4.012  49) 
+  1. 1263(0. 04838-1. 93046) 
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ACp=.84  72l  +  l2  4  8l8l-9.02l75-2  35609-  2.119  79 
ACp  =  -  .1686 

which  is  then  the  change  in  drag  coefficient  ftJi-  a  cliange  in  nominal  diameter 
from  two  feet  to  fifty  feet. 

To  estimate  the  variance  of  this  change, 


K,^  C,,  +  ’<2  ^^22  ^  ^3*^33  *^4  ^4  4  5  *^55 


+  2  K,  c,4  +  2  KiKg  c,5+2K2K3C23+2 


"2K3K3  C35+2K4KgC^5] 


.2K,K2C,2>2K,K3C,3 
+  2  1^2^  6 ‘'as*  ^  *^3*^4  ^34 


2  12  ? 

So  _0  =  006508  (6  43/' (’S)  (O  225^)M  I  lO/'/G)  (148  9(’04l 

2  ^1  ^ 

♦  (-(60  936)^(41  9722)  ♦)  3  20993)^(0  7U6  0)M  -  I  88208)^(1  56  20) 

♦  2(6  43  7  75)(  l(0726)(  2  82  30  )  •  2  ( fi  4  3  7  7  5  )(  •  I  60  9  3  6  )  (  2  I  3  5  4  ) 

♦  2(6  43775)(  3  20993  )(0  1666)*  2(6  4  3  775)(  (  88208  )( 0  3342) 

♦  2(  -  (  (O  726  )(  -  (  60936  )(  76  4  O  4  8  )  *  2  (  I  ((>  726  )(  ?  20  99  3  )(  lO  0966  ) 

♦  2  (  I  10  7  2  6  )  (  I  8820  8  )  (  I  4  96  32  )  ♦  2  (  1  6  ' )  9  3  r,  )  (  3  2  O  9  9  3  )  (  5  0  5  9  9  ) 

♦  2(  I60936)(  (8  82O8)(7  9Ol7)'2(  3  20  993)(  (eB2O8)(0  98(2)] 

-  ooesoej  094675] 

5^9  9  “  0006(6 

^2  I 

Tie-  9  ( ,ne  la  ?  1 1  ( 1 1  ■ I  a  (  i  i  in  is  ( In  -  ti  (lO'H,  I  6  w  Id  i  1 1  is  (-final  I  ■  >  ( i  O t  d  2,  and 

(In  I  y  |ii  -  I  (ill  I  lia  n  i',1  in  d  i  .1  ti  (  ( >>  1 1  i  (  i  i-nl  with  a  s  s  '  n  ia  (  cd  -  p,.  ,  t  i  nl  i  nid  i  (h  -  n  i  ■ 

limils  is  .  Md'd.  '  I  .  96  (0.  0-’.  ih7.)  Ml-  ’o.()|f<(,0.  ('Ids  i  .s  H 1 1  ( iw  n  in 

I’’  i  gn  rr-  1  ^  as  tie  s  (|na  if  s  a  lel  (  i  le  !  e  a  I  ong  w  i  I  li  lie-  o;  iiccted  cliange  in  d  la  g 

<  Of  1 1  i  (  it-n(  and  9'  pi  ?  a  i  nt  (  onlidcni  (-  limils  Iftr  any  i  liaiigf-  in  nominal  diamelei 
1 1  ( - ( w  f - f  -  n  I ' V  f  I  If  I  I  I  n  f I  M  n (  1 1 1 1  n  d  I  >  ■ ' I  If  '  -  ( 
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SECTION  m 


CONCLUSIONS 


1.  The  development  of  scaling  equations  by  statistical  techniques  is 
feasible. 

2.  Scaling  equations  were  developed  for  drag  coefficient  for  the  Cir¬ 
cular  Flat  and  Extended  Skirt  parachute  types  which  apparently  yield  reason¬ 
able  predictions  although  the  Extended  Skirt  scaling  equation  is  rather  limited 
in  range  of  application. 

3.  The  distribution  of  the  observed  data  is  probably  the  most  signifi¬ 
cant  single  factor  influencing  the  development  of  accurate  scaling  equations. 

A 

4.  A  high  multiple  correlation  (R  )for  the  regression  model  cannot  be 
relied  upon  as  a  satisfactory  condition  for  permitting  the  development  of  accu¬ 
rate  scaling  equations. 

5.  Variables,  such  as  suspension  line  length,  porosity,  etc.  ,  which 
are  known  to  affect  parachute  performance,  did  not  necessarily  appear  sta¬ 
tistically  significant  for  the  regression  analysis  due  to  the  poor  distribution 
of  available  data. 

6.  Tbe  collection  of  data  from  many  varied  sources  led  to  less  accurate 
equations  than  data  gathered  from  one  set  of  tests. 


SECTION  IV 
RECOMMENDA  TIONS 

As  noted  previously,  many  areas  in  this  investigation  have  insufficient 
or  questionable  data.  The  developed  scaling  equations,  therefore,  can  be  con¬ 
sidered  only  as  preliminary  media  in  the  development  of  basic  equations. 
Scrutiny  of  the  developed  equations  and  data  distribution  charts  and  tables  in 
Section  II  and  the  data  summary  tables  in  Appendix  II  reveals  the  lack  of  data 
in  several  areas.  There  is  practically  no  usable  information  defining  any 
parachute's  stability;  only  the  Solid  Flat  Circular  and  Extended  Skirt  types  of 
parachutes  have  usable  peak  opening  shock  and  filling  time  data.  The  Conical 
Ribbon  parachute  data  are  virtually  nonexistent.  Since  all  types  of  parachutes 
need  additional  information  to  some  degree,  the  following  recommendations  are 
made. 


1.  Future  testing  should  be  directed  toward  filling  the  existing  voids  in 
the  available  data.  The  less  the  variation  of  parameters  other  than  the  scaling 
variables,  the  more  desirable  this  testing  will  be. 
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2.  Initial  concentration  toward  defining  the  scaling  equations  for  only 
one  type  of  parachute  would  be  desirable,  since  it  is  felt  that  the  scaling  equa¬ 
tions  for  all  common  types  of  parachutes  will  be  similar. 

3.  All  future  testing  should  be  conducted  according  to  a  specific  statis¬ 
tical  design  of  experiment  in  order  to  yield  maximum  use  of  the  resultant  test 
data. 


4.  A  consistent  method  should  be  established  for  the  definition  of  para¬ 
chute  stability.  At  present,  it  is  impossible  to  correlate  the  period  and  the  am¬ 
plitude  of  parachute  oscillation  (also  gliding  and  coning)  under  free-fall  conditions 
with  wind  tunnel  tests  which  define  pitching  moment  at  angle  of  attack.  There¬ 
fore,  until  this  problem  is  solved,  pitching  moment  at  angle  of  attack  will 
probably  mean  very  little  except  in  the  general  sense  that  the  parachute  is 
stable  or  unstable  over  a  given  range.  For  an  analysis  of  parachute  scaling 
effects  for  stability,  this  consistent  definition  is  almost  mandatory. 

Looking  further  into  the  future,  the  following  recommendations  are  made; 

1.  The  range  of  applicability  of  the  scaling  equations  should  be  extended 
to  include  a  much  broader  operating  regime  (altitude  and  velocity). 

2.  The  effects  of  variables  other  than  scaling  variables —for  example, 
ribbon  spacing  —  should  be  investigated  so  that  a  scale  model  test  can  be  used 
to  evaluate  parachute  performance  without  requiring  unnecessarily  complicated 
fabrication  procedures. 

3.  The  effects  of  different  types  of  testing  — such  as  free-fall,  wind 
tunnel,  and  rocket  sled — should  be  evaluated. 
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2.  Initial  concentration  toward  defining  the  scaling  equations  for  only 
one  type  of  parachute  would  be  desirable,  since  it  is  felt  that  the  scaling  equa¬ 
tions  for  all  common  types  of  parachutes  will  be  similar. 

3.  All  future  testing  should  be  conducted  according  to  a  specific  statis¬ 
tical  design  of  experiment  in  order  to  yield  maximum  use  of  the  resultant  test 
data. 


4.  A  consistent  method  should  be  established  for  the  definition  of  para¬ 
chute  stability.  At  present,  it  is  impossible  to  correlate  the  period  and  the  am¬ 
plitude  of  parachute  oscillation  (also  gliding  and  coning)  under  free-fall  conditions 
with  wind  tunnel  tests  which  define  pitching  moment  at  angle  of  attack.  There¬ 
fore,  until  this  problem  is  solved,  pitching  moment  at  angle  of  attack  will 
probably  mean  very  little  except  in  the  general  sense  that  the  parachute  is 
stable  or  unstable  over  a  given  range.  For  an  analysis  of  parachute  scaling 
effects  for  stability,  this  consistent  definition  is  almost  mandatory. 

Looking  further  into  the  future,  the  following  recommendations  are  made: 

1.  The  range  of  applicability  of  the  scaling  equations  should  be  extended 
to  include  a  much  broader  operating  regime  (altitude  and  velocity). 

2.  The  effects  of  variables  other  than  scaling  variables — for  example, 
ribbon  spacing  —  should  be  investigated  so  that  a  scale  model  test  can  be  used 
to  evaluate  parachute  performance  without  requiring  unnecessarily  complicated 
fabrication  procedures. 

3.  The  effects  of  different  types  of  testing  — such  as  free-fall,  wind 
tunnel,  and  rocket  sled  —  should  be  evaluated. 
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APPENDIX  I 


REGRESSION  ANALYSIS  AND  COMPUTER  PROGRAM 


A.  General  Least  Squares  Multiple  Regression. 


A  general  least  squares  multiple  regression  analysis  with  k  fixed 
variates  postulates  that  some  variable  or  characteristic,  Y  ,  is  related  to 
or  depends  on  certain  other  variables  or  characteristics,  XjiXg,”*,  Xj^, 
by  means  of  the  general  linear  relation  (regression  model) 


Y=8o+^,  X,+^2 

where  X.(i  =  I,  2  ,  k  )  may  take  any  functional  form  of  the  independent  variable ; 

e.g.,  natural  log,  inverse,  etc.  Obtaining  estimates  of  the  regression  co¬ 
efficients  by  the  method  of  least  squares  consists  of  obtaining  an  expression 

Y=d;,+  b,  X,  +  X, 


such  that  the  error  sum  of  squares,  i.  e. ,  the  sum  of  the  deviations  squared 

of  the  observed  Y  from  the  predicted  or  fitted  Y,  SSE  =  ^{Y.  —  Y  be  a 

j»l  '  i 

minimum,  where  there  are  n  observations  or  data  points  (n>  k  +  |). 


In  order  that  tests  of  hypotheses  involving  the  fi's  and  the  overall 
regression  may  be  made  and  confidence  intervals  for  the  fi's  and  predicted 
Y  values  computed,  the  following  assumptions  concerning  the  measured 
variables  must  be  made: 


(1)  the  x's  are  measured  with  negligible  error  and,  therefore,  are 
regarded  as  fixed  constants  which  may  be  chosen  at  the  discretion  of  the 
experimenter: 

(ii)  for  a  given  set  of  x's  .  the  possible  Y  values  are  normally  and 
independently  distributed  with  mean,  fiy=  8q+ /3j  X|+ jSg  X2+*** +/3,,  Xj^  ,  and 
variance,  o-|  ; 

(iii)  the  variance  of  Y  for  all  sets  of  x's  is  identically  or| 

It  must  be  emphasized  that  this  postulated  relationship  between  Y  and 
the  X  's  is  not  necessarily  a  causal  one,  but  is  rather  a  relationship  which 
may  be  used  for  predictive  or  analytic  purposes. 

B.  The  Method  of  Least  Squares 

To  minimize  the  error  sum  of  squares,  the  following  procedure  is 
employed.  Write 


A  2 

vy  V* 


2 


where  the  index  of  summation  has  been  omitted  for  simplicity  of  notation. 
Partially  differentiate  SSE  with  respect  to  each  of  the  parameters  to  be 
estimated  and  equate  these  partial  derivatives  to  zero: 

■^=  2  r(V-d„-l.,X x.)(-l)  =  0 

^  =  2  Xl'^- V,  \  X.  >(-V=0 

^=2Z(y-d„-b,x-b3X,-...-b.  X.)(-X3>=0 
2  ^  • 

•  ^  • 

•  •  • 

i^  =  2Z(y-d„-b,x-b,X3-...-b.x,)(-x.)=o 

Multiply  the  first  equation  by  1/2,  perform  the  indicated  summations  term 
by  term,  and  solve  for  d^,  : 

- 

where  Y  ,  X^Xg,** are  the  arithmetic  means.  Multiply  the  remaining 

equations  by  1/2,  perform  the  indicated  summations. term  by  term,  and 
transfer  the  first  sum  in  each  equation  to  the  right  side: 

doXx,tb,Ixf  +  bjZx,x^+...+b,I  x,x,=  Zx,  y 

d;,IX2  +  b,ZX2X,  +  b2lx|+-  •  •  +  b,  I  X2X,=  Ix^y 

•  • 

•  • 

•  • 

d„Xx.+  b,Zx,x,  +  b3lx,Xj+-..-fb,ZxJ-Ix,y 


Substitute  the  expression  for  d^  in  each  of  these  equations  and  group 
similar  terms; 


b,(Zxf-x,Zx,)  +  b2(Zx,  X2-X2lx,)  +  -  •  •  +  b^(Zx,  x^-x^Ix,  )=(Zx,y-yZx,} 

b,(Zx2X,-X,Zx2}  +  b2(Zx|-X2  Z  X2  )+•  •  •+b^  (ZXgX^-X^ZXg  )  =  (Zx2Y-yZx2J 

•  • 

•  • 

*  • 

b,(Zx^x-x,Zx^j  +  b2(Zx^X2-X2Zx^)+-  .•+b^(Zx^-x^Zxj  =  (Zx^y-yZx^) 

These  equations  Ihusly  obtained  by  the  method  of  least  squares  are  commonly 
known  as  the  normal  equations. 
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Particular  note  should  be  made  of  the  quantities  in  parentheses.  These 
quantities  are  the  corrected  sums  of  squares  or  sums  of  products,  i.e.  , 

(ZXj  Xj-Xj  ZXj)=(IXiXj-  )=Z(Xj-Xj)(x^-x,)  =  (IXjX,-XjZx,), 

( ' » j  =  < .  2  ,  •  •  •  ,  k  } 

and 


(ZXj  Y- Y  ZXi)=(rx,  Y- 


(X,-X)(Y-Y),(i  =  l, £,•••,  k). 


Denoting  the  corrected  sums  of  squares  and  sums  of  products  of  the  x's  and 
Y  s  by  a|.  and  ,  respectively,  the  normal  equations  take  the  form 


<^11  ‘>.■*■“.2  ‘>2  +  ‘ 


°,k  ^  =  °,y 


"21  ‘>1  "22  ‘’2'*“ 


"2k  ^  =  "2y 


"k/^  +  «k2‘>2  +  • 


+  o 


kk 


^  = 


ky 


where,  as  indicated  above,  Ojj  =  o^j  ,  ( i,  j ,  =  / ^  2  ,  •  •  •  k )  . 

To  obtain  the  least  squares  estimates  bj ,  bg  ,  •  •  • ,  b||  of  the  regression 
coefficients  ,  0^  ,•  •  •  0 ,  it  is  required  to  solve  the  system  of  normal 
equations  for  b, ,  bg,  •  •  • ,  bj^  ,  and  then  substitute  these  values  into  the  expression 
for  dg  to  obtain  the  complete  prediction  equation  Y  • 

C.  Estimation  of  the  Regression  Coefficients  and  Tests  of  Significance 

Expressing  the  normal  equations  in  matrix  notation  where 


a.  0.-.  •  • 

•  a. 

b  ~ 

0.  If 

II  12 

(k 

1 

1  y 

A  = 

°2I  "22’  ' 

•  • 

*  "2k 

• 

,  fe  = 

*>2 

• 

,  ond  G  = 

"2y 

• 

*  • 

• 

• 

• 

•  • 

• 

• 

• 

G  a  • 

kl  k2 

’  "kk 

/k_ 

- 1 

>» 

0 

_ 1 

yields  A  b  =  G  . 

Multiplication  by  A  *  produces 

A~'  A  b  =  A“'  G 
I  b  =  a“*  G 
b  =  A  '  G  , 

where  A~'  is  the  inverse  of  the  A  matrix  and  I  denotes  the  identity  matrix. 


37 


Denoting  the  elements  of  the  inverse  of  the  A  matrix  as 


C  =  = 


C,2-- 


'Ik 


^21  ^22*  *  *  ^2k 


^kl  *^k2' 


'kk 


the  estimates  of  the  regression  coefficients  are  found  by: 


^11  “ly *'l2“2y'*'*  *  *  ^Ik  ®ky 

b2=C2,0,y +C2202y‘*'**  "^^Zk^ky 


^k”^kl®iy  ■’■^k2‘’2y‘*'“  ‘■‘■^kk^ky 


and  then 


^-‘»2  ^2 - • 

It  can  be  demonstrated^  that 


which  is  another  way  of  stating  that  the  corrected  total  sum  of  squares  for 
the  dependent  variable,  Y,  can  be  partitioned  into  a  sum  of  squares  at¬ 
tributable  to  the  regression  of  Y  on  the  x's  and  a  sum  of  squares  not  ex¬ 
plained  by  the  regression  and  which  is  presumably  due  to  the  experimental 
error,  i .  e.  , 

SST  =  SSR  +  SSE  . 


It  is,  therefore,  possible,  under  the  assumptions  of  (i),  (ii),  and  (iii),  to 
set  up  the  following  analysis  of  variance  table  to  test  the  hypothesis  per¬ 
taining  to  the  significance  of  the  overall  regression. 


1.  Anderson,  R.  L.  ,  and  Bancroft,  T.  A.,  Statistical  Theory  in  Research, 
McGraw-Hill  Book  Company,  Inc.,  1952,  p.»  171. 
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ANALYSIS  OF  VARIANCE  TABLE 


Source  of  Variation 

d.f. 

S  S. 

M.S 

F 

Regression 

k 

Ib,a.  =  SSR 

i  ' 

s|/s| 

Error 

n-k-l 

Z(Y-YJ  =SSE 
j 

Total 

n 

^(Yj-Y  )^=SST 

The  computed  F  ratio  may  be  compared  with  the  tabular  F  value  for  (k) 
and  (n-  k-l  )  degrees  of  freedom  (at  the  desired  level  of  significance)  to 
test  the  null  hypothesis  that  all  the  ^,'s  ( i  =  1 , 2 ,  •  •  • ,  k  )  are  equal  to  zero . 

Another  useful  statistic  to  be  computed  is  the  square  of  the  multiple 
correlation  coefficient, 

R^=SSR/ SST, 


which  indicates  what  proportion  of  the  variation  in  the  y's,  expressed  as  a 
sum  of  squares  of  deviations  from  the  mean,  is  explained  by  the  regression 
of  Y  on  the  x's  • 


Under  the  assumption  of  normality,  it  follows^  that  each  of  the 
bj's  ( i  =  1, 2 k)  is  normally  distributed  with  mean,  0,t  and  variance,  C|,  vf  , 
where  Cjj  is  the  i  diagonal  element  of  the  inverse  of  matrix  A  There¬ 
fore,  to  test  the  null  hypothesis  that  one  of  the  individual  regression  coef¬ 
ficients  is  zero,  h„:  j3j  =  0,  ( i  =  1 ,2  .•  •  •.  k ) ,  compute 


t  = 


and  compare  with  the  tabular  f  value  for  (n-k-l)  degrees  of  freedom  (at 
the  desired  level  of  significance).  To  test  the  null  hypothesis  that  a  par¬ 
ticular  /9j ,  ( i  =  1 , 2  ,  •  •  •  k  )  is  equal  to  some  specified  value,  compute 


t=(b. 

I 


and  again  compare  with  the  tabular  t  value  for  (n-k-l)  degrees  of  freedom 
(at  the  desired  level  of  significance).  One  other  possible  test  of  hypothesis 


1.  Anderson,  R.  L.  ,  and  Bancroft,  T.  A.,  Statistical  Theory  in  Research 
McGraw-Hill  Book  Company,  Inc.,  1952,  p.  177. 
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which  may  be  of  interest  is  the  null  hypothesis  that  two  of  the  regression  co¬ 
efficients  are  equal,  Ho:/S.=/9j(  i ,  j  =  1, 2 , •  •  •,  k  ;  i  ^  j  )  .  Testing  this  hypothesis  re¬ 
quires  the  computation  of 


t=(b, 


^(C:.+  C..-  2  Cj.)  S 


2 

E 


to  also  be  compared  with  the  tabular  t  value  for  (n-k-l)  degrees  of  free¬ 
dom  (at  the  desired  level  of  significance). 


The  fi's  are  known  as  partial  regression  coefficients  since  ySj  indi¬ 
cates  the  amount  Y  would  vary  per  unit  change  in  X.  if  the  remaining  x's 
were  held  fixed  and  bj  is  the  sample  estimate  of  this  change.  It  is  impos¬ 
sible  to  determine  the  relative  importance  of  the  various  independent  variables 
as  they  contribute  to  (or  explain)  the  variation  of  the  dependent  variable 
simply  by  examination  of  the  magnitudes  of  the  bj  s,  since  the  contribution 
to  the  overall  regression  depends  on  the  product  of  bj  times  Xj  .  However, 
the  standard  partial  regression  coefficients 

bj  =  bj  y)2  (j  =  |_2,...k) 

offer  a  proper  comparison  when  ranked  in  descending  order  of  their  absolute 
values,  since  they  are  all  expressed  in  the  same  units. 


The  upper  and  lower  confidence  bounds  for  /9j(  i  =  I  ,  2  ,•  ••,  k  )  would  be 
computed  by 

*(n-k-l)0C"/c~s| 


and  similarly  for  the  difference  between  two  regression  coefficients. 


where 
(at  the 


^i-^j,(i,j  =  l,2.---.k-,i 

L  I  \ /V  tabular 

(n-k-/  )0C 

OC  level  of  significance). 


t  value  for  (n-k  — I)  degrees  of  freedom 


To  place  confidence  bounds  around  the  average  y.Y  for  a  specified 
set  of  X  values,  ,  '  '  ’  ,  X^^^  ,  the  quantities 

^0-  \„-k-l}ccV 5'n+Z''Z  Cjj(Xoj-  X.KXq-X  )] 

i,i=l  •  -* 

are  computed,  where  n  is  the  number  of  observations  or  data  points  included 

.  ^  fc 

in  the  regression  study,  y.  =  d- +  J]  the  other  quantities  are 

o  0  *  Oj  * 

defined  as  above.  To  compute  confidence  bounds  for  an  individual  predicted 
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Y  value  given  a  specified  set  of  x‘s,  an  additional  factor  of  1  is  added  within 
the  square  brackets  of  the  above  expression  to  give 

b  J  “  ' 

The  interpretation  concerning  the  "desired  level  of  significance,  "  OC  , 
to  be  associated  with  a  test  of  hypothesis  or  with  a  confidence  interval  is  the 
following:  If  the  null  hypothesis  is  rejected  at  the  oC  level  of  significance 
(  OC  is  commonly  taken  to  be  .  05  or  .  01,  but  may  be  any  value,  depending 
on  the  problem  at  hand  and  the  use  which  is  to  be  made  of  the  results),  this 
says  that  a  value  of  the  test  statistic  (computed  from  the  sample  data)  at 
least  as  large  as  the  critical  (tabular)  value  of  the  test  statistic  would  be 
expected  to  occur  by  chance  alone  approximately  100  ‘OC  times  out  of  each 
100  repeated  samples,  assuming  the  null  hypothesis  were  true.  This  is, 
in  effect,  stating  that  the  null  hypothesis  would  be  incorrectly  rejected  100 'OC 
percent  of  the  time.  The  same  basic  idea  holds  with  respect  to  confidence 
intervals.  An  100  (1  -  oC  )  per  cent  confidence  interval  (computed  from  the 
sample  data)  for  a  given  parameter  is  interpreted  to  mean  that,  if  repeated 
samples  were  taken  and  the  corresponding  confidence  intervals  computed, 
approximately  100  (1  -  OC  )  times  out  of  each  100  repeated  samples,  the  range 
of  the  confidence  interval  would  encompass  the  parameter. 

By  application  of  the  above  computational  techniques,  it  is  possible  to 
assess  the  overall  significance  of  a  regression  model,  rank  the  independent 
variables  according  to  their  relative  contribution  to  the  total  regression,  test 
the  individual  terms  in  the  model  for  significance,  and  on  these  bases  deter¬ 
mine  if  the  model  shoxild  be  reduced.  After  a  model  involving  only  significant 
terms  has  been  obtained,  confidence  intervals  for  the  regression  coefficients 
and  predicted  Y  values  may  be  computed. 

D.  Computer  Program  for  Multiple  Regression  Analysis 

An  IBM  704  FORTRAN  II  program,  identified  as  TV-MRCA  MULTIPLE 
REGR ESSION- COMPREHENSIVE  ANALYSIS,  was  obtained  from  the  SHARE 
Library  in  Building  57,  Area  B,  Wright- Patterson  Air  Force  Base,  as  the 
most  applicable  computer  program  available.  Various  modifications  and 
additions  to  this  program  were  deemed  necessary. 

Modifications  were  required  to  render  this  program  compatible  with 
IBM  7090  FORTRAN  II  and  the  accompanying  procedures  in  effect  at  the 
ASD  computer  installation  in  Building  57,  Area  B,  Wright- Patterson  Air 
Force  Base. 

Originally,  this  routine  provided  the  least  squares  fit  of  a  multiple 
regression  model  containing  a  maximum  of  23  independent  variables  for  a 
maximum  of  500  observations  (data  points).  It  was  limited  to  the  considera¬ 
tion  of  one  dependent  variable  per  set  of  data.  The  output  included  the 
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estimates  of  the  regression  coefficients  and  their  corresponding  standard 
deviations;  the  Total,  Regression,  and  Error  sums  of  squares;  the  averages 
of  all  the  variables;  and  predictions  and  residual  errors  for  each  observa¬ 
tion  used  in  the  study. 

Two  optional  features  were  also  available.  A  list  of  synthetic  data 
points  (containing  no  value  for  the  dependent  variable)  could  be  added  to  the 
input  data  for  which  the  predicted  values  and  standard  deviations  of  the  de¬ 
pendent  variable  would  be  computed.  Any  number  or  combination  of  the 
independent  variables  could  be  excluded  from  the  regression  model  and  the 
analysis  rerun  on  this  basis  producing  the  same  output,  excepting  that  the 
computation  of  the  predicted  values  and  residual  errors  for  the  observations 
would  be  deleted. 

The  above  computer  program  was  rewritten  in  the  form  of  a  subroutine 
and  its  capacity  increased  to  provide  for  the  consideration  of  a  regression 
model  containing  a  maximum  of  49  independent  variables.  Modifications  and 
allowances  were  made  so  that  the  program  now  contains  the  following  fea¬ 
tures  and  supplements  to  its  output. 

Since  the  regression  program  is  now  in  the  form  of  a  subroutine,  a 
main  program  can  be  written  to  read  in  the  input  data,  perform  any  desired 
transformations  on  the  variables,  store  them  on  tape,  and  call  the  regres¬ 
sion  routine.  Any  of  the  variables  could  then  be  designated  the  dependent 
variable  and  a  least  squares  multiple  regression  analysis  be  run  considering 
any  number  (less  than  50)  or  combination  of  the  other  variables  as  dependent 
variables.  Thus,  the  number  of  dependent  and  independent  variables  per  set 
of  data  as  well  as  the  rerun  option  has  been  greatly  extended.  The  output  of 
the  program  now  includes  the  correlation  coefficients  between  all  pairs  of 
the  independent  variables,  the  correlation  coefficients  between  each  inde¬ 
pendent  variable  and  the  dependent  variable,  the  A  matrix  (corrected  sums 
of  squares  and  sums  of  products),  the  C  matrix  (the  inverse  of  the  A 
matrix),  an  Analysis  of  Variance  table,  the  actual  and  predicted  values  of 
the  G  vector  (the  corrected  sums  of  products  of  the  x's  and  the  Y),  the 
t  tests  for  testing  HQ*y3,=  0,  (i  =1 , 2 ,* •  • ,  k  )  ,  the  standard  partial  regression  co¬ 
efficients,  and  as  before,  the  averages  of  the  x's  and  the  y  .  the  estimated 
regression  coefficients  and  their  corresponding  standard  deviations,  and  the 

predicted  value  and  residual  error  for  each  observation  (data  point).  Values 
of  the  independent  variables  may  be  included  in  the  input  data  for  which  pre¬ 
dicted  values  of  the  dependent  variable  and  the  coi- responding  .standard  devia¬ 
tions  will  be  computed. 
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APPENDIX  II 


LITERATURE  REVIEW  AND  DATA  PRESENTATION 


A.  Literature  Review 

As  the  statistical  development  of  textile  parachute  similarity  laws  was 
based  on  the  design  and  performance  data  of  varying  sizes  of  parachutes,  all 
pertinent  literature  was  reviewed  to  gather  the  maximum  amount  of  useful 
information.  The  review  was  restricted  to  the  following  types  of  parachutes 
operating  at  indicated  airspeeds  up  to  300  knots  and  altitudes  up  to  20,000 
feet:  (1)  Circular  Flat  Solid  cloth  type  parachute  canopy,  (2)  Ebctended 
Skirt  Solid  cloth  type  parachute  canopy,  (3)  Ringslot  type  parachute  canopy, 
(4)  Ribbon  (Flat  Circular  and  Conical)  type  parachute  canopy,  and  (5)  Guide 
Surface  (ribless)  type  parachute  canopy.  Furthermore,  the  collection  of 
data  was  limited  to  the  dependent  variables  most  indicative  of  parachute  per¬ 
formance — filling  time,  opening  shock,  drag  coefficient,  and  stability — and 
those  independent  variables  considered  functionally  related  to  the  former. 

The  independent  variables  chosen  are  as  follows:  ,  Ij/D^,  Ng/D^, 

Sv/So,  X  ,W/So,  V,q,Re,M.  and  Fr. 

Only  about  a  third  of  the  more  than  150  sources  reviewed  provided 
valid  and  useful  ixiformation.  Many  bore  deceptive  titles  since  they  con¬ 
tained  insignificant  data.  Some  had  incomplete  or  unreliable  data;  for 
example,  one  source  stated  values  of  0.4  and  0.  5  for  the  drag  coefficient 
for  Extended  Skirt  parachutes.  Others  could  also  not  be  used  since  either 
most  of  the  independent  variables  were  unknown  or  none  of  the  dependent 
variables  were  given.  As  the  stability  data  in  these  sources  is  primarily 
of  an  observational  nature  and  a  statistical  analysis  requires  definitive  ex¬ 
perimental  Information,  no  useful  data  of  this  type  could  be  extracted.  All 
sources  had  to  be  considered  judiciously.  Although  similar  conditions  pre¬ 
vailed  in  some  sources,  their  data  varied  appreciably. 

Some  of  the  data  and  terms  required  modification  to  suit  the  purposes 
of  the  statistical  development.  In  some  cases  where  atmospheric  conditions 
were  not  given,  reasonable  values  were  assumed  to  calculate  dynamic  pres¬ 
sure,  Reynolds  number,  Mach  number,  etc.  The  canopy  loading  definition 
was  changed  from  gross  weight/drag  area  to  gross  weight/ surface  area  to 
provide  a  term  more  appropriate  to  the  statistical  analysis. 

B.  Data  Presentation 


All  data  are  summarized  in  Tables  6  through  10  for  reference;  these 
data  may  lead  to  additional  evaluations  in  future  studies.  These  data  are 
pertinent  to  five  types  of  parachutes  operating  at  indicated  airspeeds  up  to 
300  knots  and  altitudes  up  to  20,  000  feet. 
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Table  6.  a 

Summary  of  Data  -  Circular  Flat  Parachute 


Table  6 .  b 

Summary  of  Data  -  Circular  Flat  Parachute 


Table  6 .  c 

Summary  of  Data  -  Circular  Flat  Parachute 
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Table  6.  d 

Summary  of  Data  -  Circular  Flat  Parachute 


Table  6.  e 

Summary  of  Data  -  Circular  Flat  Parachute 
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Table  6 .  f 

Summary  of  Data  -  Circular  Flat  Parachute 
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Summary  of  Data  -  Circular  Flat  Parachute 
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Summary  of  Data  -  Circular  Flat  Parachute 
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Summary  of  Data  -  Circular  Flat  Parachute 


Table  7.  a 


Summary  of  Data  -  Extended  Skirt  (10%)  Parachute 
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Table  7.b 

Summary  of  Data  -  Extended  Skirt  (10%)  Parachute 
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Table  7.  d 


Summary  of  Data  -  Extended  Skirt  (10%)  Parachute 


Table  7.  e 


Summary  of  Data  -  Extended  Skirt  (10%)  Parachute 
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Table  7.f 


Summary  of  Data  -  Extended  Skirt  (12,  5%)  Parachute 
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Table  7.h 


Summary  of  Data  -  Extended  Skirt  (12.  5%)  Parachute 
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Table  7.i 


Summary  of  Data  -  Extended  Skirt  (14.3%)  Parachute 
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Summary  of  Data  -  Ringslot  Parachute 
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Table  8.b 


Summary  of  Data  -  Ringslot  Parachute 
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Summary  of  Data  -  Ringslot  Parachute 


Table  9.  a 

Summary  of  Data  -  Ribbon  Parachute 
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Summary  of  Data  -  Ribbon  Parachute 
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Summary  of  Data  -  Ribbon  Parachute 
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Table  9.  d 

Summary  of  Data  -  Ribbon  Parachute 
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Table  10 


Summary  of  Data  -  Ribless  Guide  Surface  Parachute 
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